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Abstract 

In this paper we study some new special functions that arise naturally 
within the framework of Hermitian Clifford analysis, which concerns 
the study of Dirac-like systems in several complex variables. In partic- 
ular we focus on Hermite polynomials, Bessel functions and general- 
ized powers. We also derive a Vekua system for solutions of Hermitian 
systems in axially symmetric domains. 
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1 Introduction 

Clifford analysis deals with partial differential operators that arise naturally 
within the context of a Clifford algebra. The operator that is studied most is 
the generalized Cauchy-Riemann operator dx + <9x, dx_ = Y^=i e jdx being 
the Dirac operator and e%, . . . , e m the generators of the Clifford algebra satis- 
fying the defining relations ejCk + ekej = —25jk- Solutions of (dx + dx)f = 
are called monogenic functions in M m+1 . One of the main basic properties is 
the Cauchy-Kowalevski extension (CK-extension) theorem which says that 
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a given real analytic function g{X_), 2L £ admits a unique monogenic 
extension f(Xg, X) is some domain of R m+1 that is given by 

f (x 0l x) = j2 L ^ L dy(x). 

k=0 

In other words, the restriction operator / — > f\x =o is injective and also 
surjective. This CK-extension leads to the construction of numerous special 
monogenic functions depending on the choice of the initial function g(X) 
and it gives rise to special classes of polynomials and other special functions 
on the way. In particular in [T7] was constructed Clifford-Bessel functions, 
monogenic generalized powers and Clifford Hermite polynomials arising from 
the CK-extension of the Gaussian distribution (the monogenic extension of 
the Gaussian distribution for m odd has been obtained in closed form in [13] 
using Fueter's theorem). The above generalized Cauchy-Riemann system is 
invariant under the orthogonal group (the spin group Spin(m) to be more 
precise) . 

There also exists a so-called Hermitian monogenic system that can be 
defined in several complex variables (see e.g. [H El El EEH EE])- It is invari- 
ant under the action of the unitary group and it contains the holomorphic 
functions in several complex variables as a subclass. To define that system 
we consider several complex variables like Zo, Z\, . . . , z n and their conjugates 
z ,zi, ... ,z n and study solutions of the system 

(f^ Z0 + fld Zl + --- + fld Zn )f(z vZj ) = (f d So +f 1 d 3ll + ---+f n d 3 Jf{z j ,z j ) = a, 

called Hermitian monogenic (/i-monogenic) functions, whereby the Clifford 
generators fj, /• satisfy the relations for a Witt basis 

fjfk = —fkfj, fj ft = — fhfji fjfk + fkfj = 

Like in the orthogonal setting one may wonder about a CK-extension the- 
orem for h- monogenic functions; it has been studied in [5]. However, this 
time the restriction operator / — > f\ Zo =o is no longer injective and it is 
also not surjective so that e.g. the Gaussian exp(— | X)j=i z jZj) admits no 
CK-extension. Yet in our paper jl] we were able to construct an analogue of 
Hermite polynomials for the unitary group (/i-monogenic setting). This leads 
to the consideration of a weaker form of the ^-monogenic system, containing 
the h- monogenic functions as a subclass, which is just big enough for the 
restriction operator / — > f\ ZQ =o to be surjective. 

The main topic in this paper is to study the CK-extension theorem for this 
sub-system of the /i-monogenic system, showing that its solutions (called h- 
submonogenic functions) are determined by their Cauchy data and to solve 
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the system explicitly for the Gaussian and for other special functions as 
well. In this way we obtain Hermite polynomials, generalized powers and 
Bessel functions for the /i-submonogenic system. We also derive a Vekua- 
type system (see e.g. [TH [18]) that describes all axially symmetric solutions 
to this system. It is the Hermitian analogue of the Vekua system for the 
so-called axial monogenic functions (see [121 EH fT7]). 

In Sections [2] and [3] we recall the details concerning the Hermitian mono- 
genic system and introduce the sub-system that will lead to a surjective CK- 
extension theorem. In Section H] we also discuss the constraints that have 
to be imposed on the Cauchy data for the CK-extension to be ^-monogenic. 
Then we derive the Hermite polynomials from the CK-extension of the Gaus- 
sian. In the next section we establish the Vekua system for axially symmetric 
solutions and apply this to the construction of generalized powers. In the last 
section we study solutions of exponential type, leading to Hermitian Bessel 
functions, which have applications in Fourier analysis. 

2 Hermitian Clifford Analysis 

In m-dimensional Euclidean space, Clifford analysis (see e.g. [31 El El E]) 
focusses on the null solutions of various partial differential operators arising 
within the Clifford algebra language, the most important one being the Dirac 
operator dx_ = Y^7=i e j^Xj , which is the Fischer dual of the Clifford vector 
variable X_ = Y^j=\ -^-j e j-> an d the null solutions of which are called monogenic 
functions. Here (ei, . . . , e m ) forms the usual orthonormal basis for the real 
vector space IR m , equipped with a bilinear form of signature (0, m) and un- 
derlying the construction of the real Clifford algebra Mo,m = ©fcLo^o,m where 
Ro m denotes the subspace of A;- vectors, spanned by the products of k differ- 
ent basis vectors. This Clifford or geometric product is governed by the well- 
known noncommutative rules e| = —1, tje^e^ej — 0, j ^ k — 1, . . . ,m. We 
refer to this setting as the orthogonal case, since the fundamental group leav- 
ing the Dirac operator dx_ invariant is the special orthogonal group SO(m), 
which is doubly covered by the Spin(m) group of the Clifford algebra Mo,m- 

When allowing for complex scalars, the same set of generators as above 
(ei, . . . , e m ), still satisfying the defining relations + e^ej = —25jk, j, k = 
1, . . . , m, may in fact also generate the complex Clifford algebra C m , where 
moreover we take the dimension to be even, say m = 2n, for intrinsic reasons, 
see e.g. [HE]. As is the complexification of the real Clifford algebra M ,2n, 
i.e. = Ro,2n © ^o,2n, any complex Clifford number A G C2 n rnay be 

written as A = a + ib, a, b e Mo,2n, leading to the definition of the Hermitian 
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conjugation: 

A^ = (a + ib)^ — a — ib, 

where the bar denotes the usual conjugation in M. 0t2n , i.e. the unique anti- 
involution for which ej = —ej, j = 1, . . . , In. This Hermitian conjugation 
then gives rise to a Hermitian inner product and its associated norm on C2 n 

given by 

(A,//) = [AV]o, |A| = >/[AtAh, 

where [ ]o denotes the scalar part. 

The above framework will be referred to as the Hermitian Clifford setting, 
as opposed to the traditional orthogonal Clifford setting. For the complex 
Clifford algebra C 2n we consider the so-called Witt basis 

fj = - ie n+j), fj = ~(ej + ie n+j ), j = 1, ... ,n. 
These Witt basis elements are isotropic 

tf = (/J) 2 = 0, j = l,...,n 
and satisfy the Grassmann identities 

fjfk + fkfj = f)fl + flf) = 0, j,k = \,...,n 
as well as the duality identities 

fifl + flfj = 5 Jk, j,k — 1, . . . ,n. 

The Clifford algebra C2 n is isomorphic to the full matrix algebra C(2 2n ) = 
End(S'), whereby the so-called spinor space S is the 2 n -dimensional minimal 
left ideal S = C2 n I, I being the primitive idempotent 

I = flfl ■ ■ ■ fnfn- 

As clearly fjl = 0, the spinor space is given by S = A n I, A n being the 
Grassmann algebra generated by the elements //,..., /t. Grassmann alge- 
bras may be seen as fermionic polynomial algebras generated by anticom- 
muting fermionic variables such as fl, . . . , /t. The corresponding fermionic 
derivatives d f u . . . , Da are clearly the elements f\, . . . , f n and the defining 
relations for <C 2 n i n terms of the Witt basis are the fermionic equivalent of 
the Weyl relations. The fermionic Euler operator Y^=i f]®f\ ls ^ ne Clifford 

number (5 = YTj=i fj fj an d ^ s eigenspaces are the spaces A e n I whereby A £ n 
is the space of fermionic polynomials of degree £ spanned by the ^-vectors 
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f\ = fat ■ ■ ■ fir The corresponding eigenvalues are given by f3A e n I = £A^I 
so that in fact j3 satisfies the characteristic equation 

P(P - 1) . . . (/3 - n) = 0. 

Using this Witt basis, the vector (X%, . . . , X 2n ) = ■ ■ ■ , x n ,yi, . . . , y n ) in 
R 2n is identified with the Clifford vector 

n n n 

K = '^2( x j e j + Uj e n+j) = 2j fj z j ~ E f^r- 

j=l j=l j=\ 

where the complex variables Zj = Xj + iyj and their complex conjugates 
~Zj = x j ~ Wji 3 = ^-1 ■ ■ ■ i n have been introduced. Defining the Hermitian 
vector variable z and its Hermitian conjugate z} by 

n n 

- = fi z ^ — = fj z ii 

j=i 3=1 

the Clifford vector X_ clearly decomposes as X_ = z — z} . This also gives rise 
to the decomposition of the traditional Dirac operator 

n n 

dx = 5>A, + en+A ) = 2 D/A " /K) = 2 (% - dj 

3=1 3=± 

in terms of the Hermitian Dirac operators 

n n 

9 * = E f} 9 * ' d ^ = = E /a 

i=i j=i 
involving the classical Cauchy-Riemann operators and their complex conju- 
gates in the complex Zj planes, i.e. d z . = \{d Xj —id yj ) and d ij = \{d Xj +id yj ), 
j = l,...,n. 

On account of the isotropy of the Witt basis elements, the Hermitian vector 
variables and Dirac operators are isotropic as well, from which it directly 
follows that the Laplacian A 2 „ = —d\ in M? n allows for the decomposition 
= 4(S £ <9 2 t + d z tdz), which is the dual expression of 

U| 2 = \?l\ 2 — + 21- 

Hermitian Clifford analysis (see e.g. [U El El HH US] ) then focusses on the null- 
solutions of both Hermitian Dirac operators obtained. Indeed, a continuously 
differentiable function g on R 2n with values in is called a /i-monogenic 
function if and only if it satisfies the system 

d z g = = d z tg, (1) 

or in other words, iff it is a simultaneous null solution of the Hermitian Dirac 
operators d z and d z \. 
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3 Introducing our system 

Let us consider the ^-monogenic system ([I]) in dimension m = 2n + 2 

(/o^o+^)/ = 

where / : fi C IR 2n+2 — >• C2 n +2 is a continuously differentiable function. If we 
multiply the first equation by /q/ (resp. fofo), we get 

Uld Z0 + ftfod*)f = (resp. / / W = 0). 

These last two equations are clearly equivalent to the first equation of ([2]). In 
a similar way, we can also show that the second equation of (J2J) is equivalent 
to 

im + /o/cfe)/ = fif d & if = o. 

In other words, / is a solution of the /i-monogenic system (J2]) if and only if 

f (/J&o + /i/oSi)/ = 

JO JO zj 

(/o^ + /o/o t ^t)/ = 
/oVo^t/ = 0. 

In this paper we shall focus on the following system of equations 



(f&o + fofodz)f = 
(/o^ () + MkO/ = 



(4) 



which is a sub-system of the h- monogenic system ([3]), solutions of which will 
hence be called /i-submonogenic functions. 

While the /i-monogenic system in C n+1 is invariant under the action of the 
unitary group U(n + 1), the /i-submonogenic system (jlj) no longer enjoys this 
invariance and it depends on the choice of the special complex direction z . 
But the system is still invariant under the action of the unitary subgroup U(n) 
of U(n + 1). Moreover, clearly the /i-submonogenic system (jlj) is equivalent 
to the inhomogeneous /i-monogenic system 



Uldz» + d*).f = fog 

(/o^o + ^t)/ = / t /i 



(5) 
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whereby g = f^d^f, h = fod z tf. This puts things back into the U(n + 1)- 
invariant setting of the inhomogeneous h- monogenic system, whereby the 
symmetry breaking comes from the specific choice of the right-hand sides 
(fog, flh) of ©. 

We will now study the Cauchy-Kowalevski extension problem for the h- 
submonogenic system. That is, given g : Q C M 2n —> C 2n+2 , find a solution 
/ of 03) such that f\ Zo =o = g- 

First note that any C 2 „+2- valued function / may be uniquely written into 
the form 

f = A + f B + f C + ftf D, 
where A, B, C, D are C 2n -valued functions. 
A direct computation then yields 

(Jo^o + fiWf = fl{d ZQ A - d z C) + flh{d Z0 B + d z (A + D)) 
(m + /o/o^t)/ = fo{d, (A + D) - d &i B) + f ft(d, f A + d- Zo C). 

Therefore the /i-submonogenic system (jlj) is equivalent to the following sys- 
tems of equations 

' d Z0 A = d z C 
d z ,A = -c\ C U 

d Zo B = -d z (A + D) 

d z ,B = d- Z0 {A + D). U 

Remark 1 Note that if moreover, we want f to be a h-monogenic function, 
then the following extra conditions should be satisfied 

d L A = d z t(A + D)= d z B = d z tC = 0. 



4 Cauchy-Kowalevski extension problem 

4.1 Power series method 

Assume 

oo oo 

A = J2( z oZo) k A k , B = J2 z o(z z ) k B k , 

k=0 k=0 

oo oo 

C = J2M z oz ) k C k , D = J2(z z ) k D k , 

k=0 k=0 



where A k , B k , C k , D k are C 2n -valued continuously differentiable functions 
denned on M. 2n . 

Remark 2 Note that f\ Zo =o = Aq + /q/o-Do- This restriction uniquely de- 
termines f. Hence, the initial conditions for the solution $8$ of the h- 
submonogenic system ^ are the functions A and D . 

Substituting the expressions ([HD in the systems §6§ and ([7]), we get the 
recurrence relations 



Ak+i 
Bk 



—d z C k . 

k+1 - 



1 



-d t tA k , k>0, (9) 



-dz(A k + D k ), A k+1 + D k+1 



k + 1 

Then, from we obtain the solution 

l -l) k 



k + l 

' d z ,B k , k>0. (10) 



k+1 



A k 
C k 



(k\y v - 



(d z d z ,) k A , 



(k + i)(k\) 2 - 

Similarly, from fflU]) we obtain that 



d z ,(d z d z ,) h A , 



k>l, 
k > 0. 



B k 
D k 



(k + l)(k\) 2 
(-l) k 



dz(d z tdz) K (A + D ), 



(k\y 



k > 0, 
k > 1. 



Remark 3 If we want the h-submonogenic CK-extension f of /| 2o =o = A + 
fofoD to be h-monogenic, the initial conditions A and D should satisfy 
the extra constraints 

d z A = d z t(A + D ) = Q. 
4.2 Clifford-Hermite polynomials 

Let us consider for example the initial condition f\ Zo =o = e ~, i.e. A = 

e~4- and D Q = 0. Note that d z A ^ 0, hence / = A + f B + $C + flfoD 
cannot be /i-monogenic. 
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Expressing A^, B^, and Df, in terms of powers of the Laplacian acting 
on the Gauss function, we obtain: 

A - (-l)*^ 1 -* )A k-i p -^ r - (~l) fc+1 4- fc fc 

r - (~ 1 ) fc+l4 '% a* „-J# 

— 1'| fc 4 1 - fe / i„|2 | ,2 



The above functions can be rewritten in terms of the Hermitian Clifford- 
Hermite polynomials introduced in jl] by means of a Rodrigues type formula 
involving both Hermitian Dirac operators: 



p G No. These polynomials can be written in terms of the Laguerre polyno- 
mials on the real line: 

HfLuJ) = C-ir 1 ^! ({n - ( l 4-) - lzz}L n p +1 



2 T 



whereby j3 = YTj=i fjfj * s the fermionic Euler operator. Hence, the functions 



Ah, Ck and take the following form: 



k (U\\2 C n 2k\Z,Z 

2- k l*l 2 / „ f\i\-\ I , ., , /'I: 



kk\ 



e 2 



(^(f)-iALri(f )),*>i 



;_l)fc+i 4 -fc 



-k-1 



B * = whm e ^ H ^ l1) = WTTy.^^ I 2 ) • k> - 
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- (¥) " ^ + (¥) (¥)) • 

where in the last line we have used the relation 

*l£2\x) = (n + a)L^(x) - nL^(x). 



4.3 Double power series method 

We now generalize the power series method of the previous subsection by 
putting 

OO oo 

A = ^ Zq z e A kti , B = Zq z e B k ^, 

k,£=0 k,£=0 
oo oo 

k,£=0 k,£=0 

with A k j, B k> £, C k j and D ki C 2n -valued continuously differentiable functions 
defined on R 2n . Note that 'f\ zo=0 = A 0fi + f B 0fi + f^C 0fi + ftf D 0>0 . 

Substitution of these expressions in the systems (JBJ) and (J7]) now leads to 
the recurrence relations 

A k +i,i = ^ ^p zPk,ti C kt e + i = — - — ^ c^t^i, (11) 

and 

B k+ \ £ = —- — — d z (A k e + Dk A, A k t + i + D k e +1 — — — — -d z tB k t. (12) 

k + 1 - « + 1 - 

Equations (|TTjl yield 

A kfi = -j-dzC k -i,o, k>l 

A k/ =(-l) k ^^-(d^) k A ^ k , l<k<£ 

A k/ = (-lY {k ^ 0,(0,t3 g )'C f fc-i-/ > o, 1 < £ < k 
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and 



C ,e = -jd^Ao^u l>\ 

C k , e = (-l) fc+l (£ " fc ^ fc)! ^t(gAt)Vi-^ 1 < k < 
C k , e =(-l) i( ^^(d^YC k ^, , l<£<k, 



while from flT2l) we obtain 



B k ,o = --j-dz(A k -i, + D k -i t0 ), k>l 



B k/ = (-l) fe ^^(9At) fc So,^, 1 < k < £ 
B k , e = (-l) m ( fc ^ 0^3^(^.1-40 + £>*-i-*,o), I < £ < k 



and 



A M + D M = (-l) fc (£ ^ fc)! a,t(aAO fc ^o,£-i- fc , 1 < A; < £ 
A M + D M = (-l)^i^ll(^ t a,)^^_ A0 + D k . lfl ), I <£ <k. 

Hence these expressions allow us to determine A kt e, B k ^, C k ^ and D k ^, when 
A 0)i , B 0;i , C kfl and D kfl are given. 

The double power series method is partitioned into three independent 
classes of solutions. 



4.3.1 Class I 



This class was considered in subsection 14.11 The functions A, B, C, D then 
take the form (El). 



4.3.2 Class II 

In this case we put 

oo oo 

A = 4 J^^oM, B = 4 +1 J2(z z ) k B k , 

k=0 k=0 i q\ 

oo oo \ 16 ) 

C = 4~ l Y,( z ^o) k C k , D = z s ^2(z z ) k D k , 

k=0 k=0 
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with seN and A k , B k , C k , D k C 2n - valued continuously differential) le func- 
tions defined on M? n . 

Remark 4 In case of s > 1 we have that f\ ZQ =o = 0, while in case of s = 1 
we obtain f\ Zo =o = foCo- However the restriction f\ Zo =o does not determine 
the function f uniquely. 

Plugging expressions ( TT3j) in the systems of equations §6§ and (J7j), we obtain 
the following recurrence relations 

A k = -^—d z C k , C k+1 = --^—d zt A k , k>0 (14) 
s + k ~ k + 1 ~ 

and 

B k = -J— r- d z (A k + D k ), A k+1 + D k+1 = ——d z1 B k , k > 0. (15) 

s + 1 + k ~ k+1 - 

Next, ([HD leads to 

A k = ("i) fc ^|y^t^)^o, k > o 

while using (fT5j) we get 

D * = (-1) ^(^)1 (s{d,,d,) k D - d^dtfCa) , k > 1. 

Hence, the above equations allow us to determine all A k , B k , C k and D k 
given the starting values C and .D . 

If we take for example as initial conditions the Gauss function, i.e. Co = 

_|z|2 

D = e 2 , the solutions can again be expressed in terms of the Hermi- 
tian Clifford-Hermite polynomials as before or, alternatively, in terms of the 
Laguerre polynomials: 
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j2 

k\(s + l + k) ° 

s! Ul 2 i / 1-2 



(s + 1 + fc)! — 



2fc + ]A^' 2 



-^-*{f((-.-^(f)-^ /U12 



> i. 



4.3.3 Class III 

Here we look for solutions of the systems ([6]) and (J?D of the form 

oo oo 

A = r o ^(z -z ) k A k} B = z s Q - 1 J2( z oZo) k B kl 

k=0 k=0 

oo oo 

C = Zq +1 y^(^o^o) fc Cfc, D = Zq }](z z ) k D k , 

k=0 k=0 



(16) 



with again s£ff and A k , B k , C k , D k C2 n - valued continuously differentiate le 
functions defined on M. 2n . 



Remark 5 We have that 

/I ^0=0 



ifs>l, 
f B ifs = l. 



However, again this restriction does not determine the function f uniquely. 

Substituting the expressions ffTBT) in the systems (J5J and ([7D, we arrive at the 
following recurrence relations: 

A k+1 = -^—d,C k , C k = —J— -rdjA k , k>0 (17) 
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and 

B k+1 = - T —d,{A k + D k ), A h + D k = -J—SgtB*, A;>0. (18) 
From (j!7p we obtain expressions for and in terms of A$\ 

Ct = ( - 1)t+ ' fe!(s+ t + 1 ) ! ^ (a A ) ^ 

while fll8p yields expressions for B^ and in terms of the starting values 
Aq and Bq\ 

Let us again illustrate the above with the Gauss function as initial conditions, 

\z\ 2 

i.e. A = B = e 2 . Like before, the solutions can be expressed in terms of 
the Hermitian Clifford-Hermite polynomials, or, alternatively, the Laguerre 
polynomials: 



Sfc - ( " 1)4 ^!( s + fc-l)! e 2 ^,U^) 

A;( S + A;-l)! e ^ 2 J 2" V 2 



^•^^ ; feu f¥l - f¥ ) ) , * > 1 



2" fc "S S ^-— Vy e-^ zU k l ( % ) , fc > 

(s + k + l)l ~ V 2 
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k A -k ( s ~ X ) ! „-JC (itO) 



l 2fc+l 



(s + ife)! [- fc V 2 ; V V 2 

1 + (\Z\ 



))}' *^ 



2 

D n = -e-^|-, + l 

5 Axial-type solutions 

5.1 The Hermitian Vekua system 

Here, we will look for special solutions of the systems (0) and (J2J), which we 
assume to be of the form 

A = a\ + Ao2, .B = Zoz}b, 

C = zqzc, D = di + zjzd 2 , 

where ai, a 2 , b, c, d±, d 2 are continuously differentiable functions depending 
on the two variables 

(y ,u) = (\z \ 2 , \z\ 2 ), 

and taking values in the real algebra generated by (5 = YTj=\ fjfji which, as 
we know, is the quotient of the free ring generated by /3 with its 2-sided ideal 
generated by the polynomial fi(/3— 1) . . . (/?— n). One may also look for spinor 
valued solutions in K l n I thus allowing (3 to be replaced by its eigenvalue I. 
It is easily seen that 

d ZQ A = z (<9^ai + z}zd U0 a 2 ) 
d & C = z ([d,z)c + J2 J = *° O + 

&z Q C = z(c + v d uo c) . 
Using the identity 

zf3 — f3z = z, 

we also obtain that 

n 

dz?A — zd v a\ + fj (fjS. a 2 + Zjz}zd v a2j = zd v a x + (n — /3)za 2 + zz}zd u a 2 
i=i 

= z {d v a\ + (n + 1 — /3)a 2 + vd v a 2 ) . 
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Therefore, system ([HD takes the form 

d UQ ai = /3 c 

d vo a 2 -d u c = (19) 
d u ai + vd v a 2 + vod Uo c = (/3 — n — l)a 2 — a 

Making similar calculations, we may check that system ([7]) can be rewritten 

as 

d v (ai + di) + u d Vo b = f3(a 2 + d 2 ) - b 
d^iat + di) -vd v b = {n-/3)b (20) 
d Uo (a 2 + d 2 ) + d v b = 0. 



Remark 6 Note that systems nW\) and §20y are the Hermitian equivalent of 
the Vekua system considered in JX^, ITSj [77] / that describes axially symmetric 
monogenic functions in the orthogonal setting. 

In order to solve the Vekua systems ([191) and (|20|) . we shall make use of the 
power series method by writing 

oo oo 

k=0 k=0 

oo oo 

fc=0 fc=0 

with j = 1,2. The above systems may now be rewritten in the form 

a k ,i = yc^i, k > 1 (21) 



k 



a k>2 = ^=i, k > 1 (22) 

k 

(k + l)c k + af kl + ua' k2 = (P-n-l)a k2 , k > (23) 



(k + l)b k + a' kA + d' kjl = f3(a k>2 + d k , 2 ), k > (24) 

1 
k 



a k ,i + 4,i = z((n- /3)6 fc _i + vb' k _^) , k>l (25) 



a k , 2 + d k>2 = -^, k>l. (26) 
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Substituting (ED, fl22J) in (1251), we get 

fc(Jfe + 1) v ; ~~ (27) 

Co = {P - n - l)a ,2 - « ,i _ w d,2- 

In a similar way, substituting f|25|) . f }26|) in f|24|) . we obtain that 

= "MITT) «" + + "*->) ' * S1 (28) 

6 = /3(ao,2 + ^0,2) — a. 0)1 — d 01 . 

Finally, from (p]) and using Q2B, (|22J we get 

4.1 = r O - 0)&fc-i + i/^-i " /9cfc-i) , fc > 1 (29) 

4.2 = ~(65b_i + 4-i). fc > L (3°) 

Hence the explicit solutions of the Vekua systems (TT9l) . (120]) may be obtained 
via the recurrence relations Q2ZJ), ([28]), (EI]), ((22]), (]29]) and (gQJ with the 
initial conditions 

%(o,^), 4,(0,1/), j = 1,2. 



5.2 Hermit ian generalized powers 

Let us illustrate the above ideas with an example that relates to the definition 
of generalized powers in the Hermitian setting. For the orthogonal setting 
we refer the reader to [17J). Suppose that 

«o,i(^) = «i^ s , a 0t2 (v) = a 2 v s ~ 1 , 

da,i(v) = 5xu s , d 0t2 (u) = <5 2 z/ s_1 , 

with s,a>j,8j G M, j = 1,2. Using the previous recurrence relations we easily 
obtain 

b k (v) = a(k)u s ~ k ~\ k > 0, 
Ck (v) = S(k)i/ S ' k ~ 1 , k > 0, 

where 

:-i) fc nti(^ - q nti(n + * - *) (/% 2 + 52) - («i + w 



o;(fc) 



(fc + l)(fc!) 2 
17 



_ (~l) fc Uti(s - I) nL(n + g - g ((g - n - s) Q2 - ai s) 

{> (k + i)(k\y 

Hence 

. . p8(k-l) s _ k 
a k ,i{v) = tL ~^ V k , 

_ {s- k)5{k - 1) k _ x 
ak,2\v) — v , 

(7i + s-fc-/3)a(fc-l)-^(fc-l) 
dkM = "^X^ " X ) + ^ - !)) Z/S ^ 1 ' 

with k > 1. 

Remark 7 zs worth pointing out that if s G N £/ien a(fc) = = 
/or k > s. Therefore, for the case s G N £/ie solutions A, B, C , D are 
homogeneous polynomials of degree 2s in M 2n+2 . 



5.3 Two more cases 

We may consider two more cases of axial-type solutions of the systems (jH]) 
and (j7j). Indeed, if we assume that 

A = zg(ai + -B = -?o + V fe ' 

with s G N, then we obtain the following systems of equations 

vod Vo a 2 - d u c = -sa 2 
d v a\ + vd v a 2 + d Vo c = (f3 - n - l)a 2 

d u (ai + d\) + V(jd VQ b = j3(a 2 + d 2 ) - (s + l)b 
d vo (ai + di) - vd v b — (n — /3)b 
d V0 (a 2 + d 2 ) + d v b = 

which can be solved using the power series method as we did in subsection 
15.11 The explicit solutions are given as follows 

c * = ~ k{k + \ _ 1} ((" + + <-x) , fc>l 
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ak,i = j^— c k , k>0 
k + s 



>>*=- k(k+ \ +1) ((n+m- 1 +<- 1 ), *>1 

fo o = ^-^-(^0,2-^0,1) 
4,i = ^ ((n - P)b k -i + ^Li) - £+7 °k, k>l 



k k + s' 
with the initial conditions 



4,2 = - fc > 1 



c(0,i/), 4(0, z/), J = 1,2. 

Finally, let us assume that 

A = Zq(cli + z}za 2 ), B = z^z^b, 

C = z s +1 zc, D = z s (di + z}zd 2 ), 

with s G N. This yields the systems 

d„ ai = f3c 
d U0 a 2 - d v c = 
d v a x + vd v a 2 + v Q d v ^c — (/3 — n — l)a 2 — (s + l)c 

d v (ai + 4) + d VQ b = /3(a 2 + 4) 
vodvo ( a i + 4) - vd v b = —s(a 1 + 4) + (n — f3)b 
vod vo (a 2 + 4) + d u b = -s(a 2 + 4), 
whose solutions may be obtained via the recurrence relations 

Ct = - k(k + \+i) + <-.)■ *si 



c 



1 -((/?- n - l)oo,2 - a ,i - ^0,2) 



s + 



Ofe,l = ^ Cjfc_i, fc > 1 

c' 

«fc,2 = > 1 
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k(k + s - 



((n+l)b' k _ 1 + ubU), k>l 



6' 

4,2 = -77 a* 2, fc > 

with the initial conditions 

0,(0,1/), 6(0,1/), j = 1,2. 



6 Exponential-type solutions and Hermitian 
Bessel functions 

Here we shall seek special solutions of the systems and ([7j) of the form 

A = e Xzo+ ^° (ai + z}za 2 ) , B = e Xz ° + ^°z}b, 

C = e Xzo+ » z °zc, D = e Xzo+ ^° (4 + z}zd 2 ), 

X, fj, G R \ {0} and where a±, a 2 , b, c, 4, 4 are continuously differentiable 
functions depending on the variable v — |z| 2 taking values in the real algebra 
generated by j3. 

In what follows J a , I a stand, respectively, for the Bessel function of the 
first kind and the modified Bessel function (see e.g. [ID]), given by 

k / + \ 2k + a 



j it) = y (t) 



oo j , x 2k+a 

A;!T(A; + a + 1) \2 



Proceeding in the same spirit as in the previous section we can check that 
systems ([6]) and (J7|) now take the form 

ai = -c 
d 

a2= x 

^ a[ + va! 2 + /ic = (ft — n — l)a 2 
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a[ + d[ + Xb = (3(a 2 + d 2 ) 

ai + dx = - ((n - (3)b + vb') 
/i 

b ' 

a 2 + d 2 = . 

H 

From these systems of equations we can easily deduce that b and c satisfy 
the same ordinary differential equation, namely 

vb" + (n + l)b' + X^ib = 
vc" + (n + l)c' + A/xc = 0. 



It has solutions of the form 



b(u) = I aiU ( 2 v(V^)_ if V 
^ a,\V~^I n (2y/—\/jv) if \ji 

c(v) = 



>0, 
<0, 



^J n {2y/XfjjPj if A/i > 0, 
a 2 v~2l n if A/i < 0, 

where a±, a 2 are arbitrary real constants. The other functions can be easily 
obtained from the above equations. Indeed, 



a 2 {v) 



y/Xu _n+l 



a 2 



+ i if A/x > 0, 

V^-.-=*!r ( 2v ^V;) if A/i<0, 



A 



Z/ 2 J, 



n+1 



-V 73 ^ + ^"^4+1 (2 V /Z V^) if A/i < 0, 



d 1 (u) = -(n - P)b(v) - L(v) - v(a 2 (v) + d 2 (v)). 

/i A 
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